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ABSTRACT 


Historically, studies on decorrelation scales have been conducted in the deep 
ocean waters. As the Navy shifts its interest toward the less understood shallow 
water regions, decorrelation scales need to be computed in order to use formerly 
deep water models such as the Optimum Thermal Interpolation System (OTIS) for 
shallow water regions such as the Yellow Sea. 

A data set containing over 35,000 temperature profiles from 1929 to 1991 
was obtained from the Naval Oceanographic Office’s MOODS data set. The winter 
and summer seasons provide realistic results. Winter has the smallest decorrelation 
scales of all the seasons, approximately 15 days and 165 km. Summer shows that 
there are different decorrelation scales between the surface and at depth. The surface 
has scales of 12.3 days and 251 km while at depth the scales are approximately 16.5 
days and 163 km. An observational sampling network design is suggested for future 
sampling of the region. Spring and fall provide mixed results which may be due to 
the irregularities in time and space of the data set or to the very complex forcing 
mechanisms found in the region. Overall, this study gives a ground work for better 
refinement of decorrelation scales and thus, the ability to assess the conversion of 


deep water models to shallow water regions. 
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I. INTRODUCTION 


Various ocean systems such as fronts, eddies and water masses have different 
temporal and spatial correlation scales (Figure 1.1). These scales feature a system’s 
life span and spatial extent both horizontally and vertically. Horizontal scales are 
often divided into north-south and east-west length scales and are represented by a 
variety of parameters. For example, Provost and Le Traon (1993) attempted to 
describe the Brazil-Malvinas Current confluence confines through temporal and 
spatial scales using variable altimetric measurements. 

Horizontal scales for temperature fields in the deep ocean basins have been 
extensively studied. Each ocean has its own general scales with variations in specific 
regions. White et al. (1982) identified spatial correlation scales in the western 
Pacific of about 600 km in the tropics (south of 17.5°N) and 300 km in the subtropics 
(north of 17.5°N) while Ozsoy et al. (1989) found the spatial scale to be 200-250 km 
in the Levantine Basin of the Mediterranean Sea. Historically, most of the studies on 
correlation scales have been conducted in deep ocean waters; few have been 
conducted on the continental shelf. There are several reasons for this occurrence 
with the main one being that deep water physics are better understood. Another, yet 
not an obvious reason, is that in past years the US Navy has put its emphasis on deep 
water operations. Therefore, most of the Navy’s research and modeling assets have 
been directed toward deep water applications. This poses many challenging 
problems as the Navy shifts its interests toward the less understood shallow waters 
regions. Naval models such as the Optimum Thermal Interpolation System (OTIS) 
that have been designed for deep water usage need to be modified for shallow water 
use. 

Deep ocean temperature profiles tend to follow the pattern of a mixed layer, 


a thermocline and a deeper layer slowly decreasing in temperature with depth. 
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Figure 1.1 Approximate Ranges of Horizontal and Temporal Scales of 
Selected Ocean Variabilities (from Mooers and Robinson, 1981). 


Shallow water does not consistently mimic this pattern. It may range from the 
classical profile or be completely isothermal. The shallow environment influences 
the temperature distribution with some parameters like bathymetry, river run-off and 
tides, often omitted in deep water, becoming increasingly important in shallow 
water regimes. Spatial scales are strongly controlled by bathymetry. In shallow 
water, besides scales associated with seasonal changes, some other time scales also 
become important. The forcing of internal waves (minutes), tides (hours-day), 
meteorological forcing (hours-weeks), freshwater runoff (days-months), and 
interactions with slope or deep water (days-weeks) can have a great impact on 
temporal variability (Haeger, 1991). 

The Yellow Sea was selected for temperature horizontal correlation analysis 
based on two main considerations. The first is that the Yellow Sea demonstrates 
shallow water characteristics for which a data set of 50,000 profiles exists. The 
second consideration is from a military standpoint. The region is a potential “hot 
spot’ and needs special study. Modification of the OTIS model for the Yellow Sea 
can be used for other marginal seas. To do so, there must be an accurate climatology. 
Since this is a new region of emphasis, as stated earlier, a climatology has not been 
built from the Naval resources. Thus, there are two tasks in this study, (1) 
constructing a seasonal climatology from the data base and (2) computing horizontal 
and temporal scales and their vertical variation. This work will provide useful 
information for the OTIS model, since it needs temporal and spatial scales prior to 
running the model. After the climatology has been established, temperature 
anomalies are computed by subtracting the climatology from the observations. The 
anomalies are then paired and placed into bins depending on their separation in time 
and space. A correlation coefficient is then calculated from all the anomaly pairs in 
each bin. An autocorrelation coefficient matrix will be constructed for the 


combination of different lags. This matrix represents the temporal and spatial 


variations of the autocorrelation function which is fitted by the Gaussian function 
used in the OTIS model. Finally, the temporal and spatial correlation scales are 
obtained. 

To investigate the vertical variation of the horizontal and temporal scales 
several schemes have been used in deep water analysis. Sprintall and Meyer (1991) 
used the sea surface temperate (SST) and the depth of the 20°C isotherm for the 
Eastern Pacific and White et al. (1982) used 100, 200, and 300 m depths for the 
Northwestern Pacific. Neither scheme is practical for the Yellow Sea region since 
the 20°C isotherm is not always present and the water depth is less than 150 m 
throughout. Since shallow water is strongly influenced by bathymetry, it was 
decided that the seasonal correlations scales for the temperature fields should be 
computed for the surface and depths representing 50% and 80% of the water 


column. 


Il. REGIONAL CHARACTERISTICS 


A. GENERAL 

The Yellow Sea is the middle section of a continental sea system (Figure 2.1). 
To the north is the Bohai Sea and to the south the East China Sea and the Western 
Pacific beyond. It lies between mainland China and the Korean peninsula and covers 
an area of approximately 404,000 km. There are four major rivers that flow into it. 
The Yangtse River located to the southwest, the Yellow and Liao Rivers to the north 
and the Han River to the east. 

The Yellow Sea is considered a shallow basin sea. The bottom slopes from 
20 m in the north to its southern boundary, the shelf break, at 150 m. The Bohai Sea 
region ranges in depth of 20 to 50 m. A trough oriented north to northwest and 
paralleling the Korean coastline occupies the center of the sea. It has a maximum 
depth of 90 m. Within 50 km of the coastline the average water depth is 20 m. 
Moving toward the shelf break the depth slowly increases. The shelf break runs in a 
roughly northeast-southwest orientation and quickly drops from 150 to 500 m in a 
few kilometers. For the remainder of this study the Yellow Sea will include the 
Bohai Sea and that portion of the East China Sea that lies on the continental shelf 
unless specifically stated. 


B. METEROLOGY 

Nester (1977) compiled an extensive report that is the foundation of this 
meteorological overview with wind speeds supplemented by the US Navy Regional 
Climatic Study (Elms, 1989). The Yellow Sea lies within the Asian monsoon 
circulation pattern and experiences hot, short, humid summers and long, cold, dry 
winters. Winter is dominated by an intense cold high pressure system called the 
Siberian high. Its greatest intensity occurs in January, pushing cold, dry air from the 
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Figure 2.1 Bathymetry of the Bohai, Yellow and East China Seas. Depths are 
in 10 m Increments. Bathymetry Data are from the NGDC and 
NAVOCEANO World Wide DBDS5 Bathymetry Data Bases. 


interior of the continent along the east coast of the Asian mainland. The clockwise 
circulation of the Siberian high produces northerly or northwesterly winds referred 
to as the northwest monsoon. This is the predominant feature of the winter weather 
even though small iow pressure areas known as Shanghai lows form over the Yellow 
Sea and travel over the Korean peninsula. Figure 2.2 shows the mean pressure and 
wind systems for February. Surface winds are normally from the north-northwest 


with an average speed of 8 m/s. 
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Figure 2.2 February Mean Pressure and Surface Wind Pattern During the 
Winter Monsoon (from Langill, 1976). 


March marks the beginning of the transition from winter to spring. The 
weakening of the Siberian high permits the northerly migration of the polar front. 
Tropical maritime air begins to push northward displacing the cold, heavy Siberian 
air causing “migratory highs” that may in turn be displaced by low pressure areas 
from the south. Northerly winds are no longer dominant and by the end of the season 
winds are from the south averaging 4.6 m/s. 

Summer begins in June and is also referred to as the southwest monsoon 


season. The polar front has migrated to its farthest northern extent. It is a boundary 


zone between continental! air masses from Asia and the warm motst air from the 
tropical oceans and serves as a buffer zone between the two air masses. Figure 2.3 


shows the mean pressure and wind systems for the onset of summer. As summer 
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Figure 2.3 June Mean Pressure and Surface Wind Pattern During the 
Summer Monsoon (from Langill, 1976). 


progresses, the Bonin high pressure cell expands westward into the East China Sea 
while the Manchurian low, caused by land mass heating, settles over northern China, 
just north of the Yellow Sea. This is the season of predominantly southward winds 
and heavy rains. Although the three months of summer are ripe for typhoons, the end 
of July and the beginning of September is the period of their greatest probability. 
The East China and central Yellow Seas have a 74% and 32% chance, respectively, 
of a tropical cyclone occurring in the region (Elms, 1989). The polar front ts unstable 
and 1s conducive to heavy rainfall. The season starts with prevailing winds from the 
south at 5.1 m/s but shift to the north-northeast by September. 

September marks the fall transition season with a rapid southward movement 


of the polar front accompanied by a wide band of rain. By the end of the month the 


southerly warm, moist, tropical air of the summer monsoon has been replaced by a 
cooler, drier air from the continent. The Siberian high begins to form as the Asian 
continent cools causing cool air masses to periodically break away and be carried 
southward. Winds are predominately from the north at 8.2 m/s. Table 2.1 


summarizes the seasonal cycles of weather phenomena. 


C. OCEANOGRAPHY 

The seasonal changes of the Yellow Sea are strongly influenced by the 
Southeast Asian monsoon. Nester (1977) discusses this monsoonal influence on 
SST, mixed layer depth (MLD) and current speeds and 1s the basis for this section. 
In winter, when the monsoon winds are at their strongest, the Yellow Sea Current is 
also at its strongest, flowing southward averaging 51 cm/s to the end of the Korean 
Peninsula (Figure 2.4a). Hsueh (1988) found that high northerly wind events during 
winter create a mean sea level gradient that forces bottom water to flow northward 
in spurts along the Yellow Sea trough. Sea surface temperatures in the Bohai Sea 
cool to the freezing point while in the Yellow Sea it is 6°C at its northern extent and 
10°C at the shelf break. The MLD is at its deepest during winter owing to both 
convection and wind mixing by the strong monsoon winds. 

March marks the weakening of the northerly winds which in turn leads to the 
weakening of the Yellow Sea Current. The current is variable in direction and speed 
until the end of the season, when it starts flowing northward. The sea surface slope 
is reestablished and a warming trend of the Yellow Sea Current begins (Hsueh, 
1988). The SST’s quickly transition from winter to spring temperatures with an 
average increase of 10°C. The Bohai Sea starts at 6°C and by the end of the season 
it is up to 16°C which is similar to the rest of the Yellow Sea. The MLD is in 
transition so both shallow and deep layers can be found. This is the season of least 


wave height in the Yellow Sea. 
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Table 2.1: Climatic Calendar for South Korea (from Langill, 1976). 
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By the time July arrives the Yellow Sea Current has developed into a closed 


cyclonic circulation (Figure 2.4b) which causes an upwelling effect in the central 
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Figure 2.4 Representation of Surface Currents Off the Korean Peninsula in (a) 
Winter and (b) Surnmer (from Bartz, 1972). 


region. This circulation occurs above the thermocline, which is about 13 m in 
August, and may reach a speed of 15.5 cm/s. Below this circulation is a cold water 
mass, commonly referred to as the ‘Cold Water Mass of the Yellow Sea’ (CWYS), 
that remains unchanged and motionless throughout the summer (Li and Yuan, 
1992). The SST is nearly uniform throughout with an average temperature of 23°C. 
The MLD is at its shallowest and maybe negligible in the coastal and the Bohai Sea 
regions due to surface heating. The water is strongly stratified. 

October is the beginning of the transition back to winter conditions. The 
southerly winds have weakened letting the sea surface slope reestablish itself again. 
The Yellow Sea warm current becomes weak and variable. The SST steadily 


decreases 20°C from the beginning of the season to the end. The MLD begins to 


I} 


deepen as stronger winds and convective mixing occur, yet the water mass is still 
Stratified. 

Liu et al. (1992) presents a loose interpretation of the temporal and spatial 
distribution of water masses through temperature and salinity analyses. The Yellow 
Sea consists of five main water masses: the coastal water, the Yellow Sea surface 
water, the Bohai- Yellow Sea mixed water, the Yellow Sea bottom cold water and 
the Kuroshio water. The coastal water is influenced by the amount of river runoff. 
In the spring and summer, when run-off is the greatest, low salinities (28.4-30.5) and 
high temperatures (27.5°C) identify this water with the reverse true in the fall and 
winter. The Yellow Sea surface water encompasses the middle of the sea. It 1s the 
dominant water mass in the summer due to the intense surface heating. During the 
fall and winter months the water cools and becomes one water mass with the bottom 
water. The Bohai- Yellow Sea mixed water is in the northern region where the two 
seas connect. This water mass is a mixture of the coastal waters and the Yellow Sea 
surface water. Its greatest extent is in fall and winter when it reaches the Korean 
Peninsula. By summer the Yellow Sea surface water has penetrated into the Bohai 
Sea. The Yellow Sea bottom cold water or the CWYS is present in the spring and 
summer. Its mean temperature is 5 to 12°C and the mean salinity 31.8-33.8. As 
stated earlier, with the approach of winter the CWYS merges with the surface to 
form one water mass named the Yellow Sea surface water. The Kuroshio water is 
located along the shelf break and is characterized as high salinity and high 
temperature. The mean temperature is 18-29°C and the mean salinity is 34.2-34.9. 
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Iii. OPTIMUM THERMAL INTERPOLATION SYSTEM 


The US Navy requires and conducts daily analyses of the ocean thermal 
structure for a wide variety of tasks. These analyses are made for both global and 
regional areas, with resolutions as small as 20 km in some regions. The Navy’s 
Optimum Thermal Interpolation System (OTIS) combines real-time data, 
climatology and predictions from ocean mixed-layer models to represent an 
accurate picture of the ocean thermal structure on global and regional scales 
(Phoebus,1988; Clancy et al., 1989). OTIS can be executed in either “shallow” or 
“deep” analysis modes. This statement is somewhat deceptive since the term 
shallow in this case refers to a depth of 400 m. The analysis sequence for OTIS is 
the calculation of (1) SST, (2) the mixed layer depth, and the (3) temperature below 
the mixed layer depth. In truly shallow waters such as the Yellow Sea this analysis 
sequence may be inappropriate since isothermal profiles are possible. 

OTIS is designed to receive its climatological data from the Generalized 
Digital Environmental Model (GDEM) which synthesizes a climatology from a 
given data base. Adequate climatologies for the world’s marginal seas are not 
prevalent but are being rapidly developed by the Naval Oceanographic Office. Even 
though OTIS is not a shallow model in the literal sense, it is hoped that with the 
proper parameters OTIS can be modified for use in marginal seas. 

OTIS is based on the Optimum Interpolation (OID) data assimilation technique 
of Gandin (1965). OI combines a climatological temperature at a gridpoint plus a 
weighted combination of observed and model predicted anomalies. Thus the 


analyzed SST at the k th gridpoint is 


N 
T, = 1, + >) Oyj (T; -~ 7) +B, (4-7) G.1) 


i= 1 
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where 7;,° is the climatological surface temperature at the gridpoint, 7,? the predicted 
temperature at the gridpoint produced by a 24-hour forecast, T° the observed 
temperature at location i, 7; the climatological temperature at location i valid at the 
time of observation i, a, the weight applied to the observational correction, B, the 
weight applied to the prediction correction, and N the number of observations. With 


the minimal mean square error, these weights satisfy the following equation: 


N 
ll (ge 5,45) Ops + (M1, +3;,44) By Exe) 
fl 

where 1, is the space-time autocorrelation function between climatological 
anomalies at location i and location j, A? is the noise-to-signal ratio for the 
observation, i,’ is the noise-to-signal ratio for the prediction and 6, 1s the Kronecker 
delta function. The autocorrelation function 1s represented as 

N; j 
where Ax;, Ay, and Az, are the east-west, north-south and time separations between 


ee 2 2 re 7 
= exp(-A’ Ax’ ,-B’,Ay p-C At ;) (3.3) 


locations. The quantities A,, B, and C, are the inverses of the east-west, north-south 
and time decorrelation scales at different depth levels, k. These parameters are user- 
specified within OTIS and must reflect the oceans’s resolvable variability about 
climatology. The better the correlation scales are calculated, the higher the accuracy 
the thermal field will be. In order to implement OTIS to the Yellow Sea region 


decorrelation scales, Va By! and C, ; , must be determined. 
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IV. METHODS OF ANALYSIS 


A. MASTER OCEANOGRAPHIC OBSERVATION DATA SET (MOODS) 
OF THE YELLOW SEA 


1. General 

The oceanographic data used for this study were taken from the Master 
Oceanographic Observation Data Set (MOODS), the U.S. Navy’s environmental 
data base. It is a collection of observed vertical profiles of temperature, salinity and 
computed vertical profiles of sound speed. MOODS consists of over five million 
world-wide profiles of oceanographic data and is considered to be the most 
comprehensive collection of temperature and salinity data known to exist. The data 
base has been compiled from military, academic and commercial sources, both 
foreign and domestic. It has a classified and unclassified version. The data were 
derived from various instrumentation types ranging from the Nansen bottles of the 
1920’s to the present day Conductivity-Temperature-Depth (CTD) recorder. The 
MOODS data are not in a gridded form; they are irregularly distributed in space and 
time. 

The original MOODS data base used for this study consisted of 
approximately 50,000 temperature profiles located between 117-128 °E longitude 
and 25-41°N latitude (Figure 4.1a) which encompasses the Yellow Sea and the 
northern East China Sea. The data spans the period from 1929 to 1991. The data 
were screened to restrict observations only to those temperature profiles that were 
located on the continental shelf. This reduced the data set to 35,658 profiles (Figure 
4.1b). Since horizontal scales have seasonal variations, the data base was binned into 
four seasons. The seasons were defined according to the convention of the Naval 
Oceanographic Office for the Yellow Sea (Haeger, personal communication) 


January, February and March constituted winter; April, May and June, spring; July, 


I5 


August and September, summer; October, November and December fall. Table 4.1 


shows the partition of the data; winter has the least profiles and summer the most. 


(a) —— (b) 
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Figure 4.1 (a) Distribution of Entire MOODS Temperature Profile Data Set, (b) 
Profile Distribution After Removal of Data Deeper than Restricted to the Shelf 
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TABLE 4.1: Number of Temperature Profiles in the Data Base for Each Season. 





The number of temperature profiles from year to year fluctuates greatly. The 
period of 1963 to 1977 is found to have the greatest concentration of profiles 


averaging 315 per season. 
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2. Limitations 

The main limitation of the MOODS data is its irregular distribution in time 
and space. Certain areas 4nd events are over sampled while other areas lack enough 
observations to gain any meaningful insight. Vertical resolution and data quality are 
also highly variable depending much on instrument type and sampler expertise. 

The sparsity of profiles have a direct impact on the type of analysis conducted 
on the data. A prime example of this is the Bohai Sea region in winter. Access for 
data gathering in this region has been limited. In fact, most observations in the area 
have been conducted prior to the 1960’s. Figure 4.2a depicts unclassified 
observation locations for the winter seasons of 1963 to 1977. Another data poor area 
is the eastern coastline of mainland China (Figure 4.2b). The spring of 1969 had the 
greatest number of profiles for a given season yet the entire Chinese coastline lacks 
a single recorded hydrographic cast. An analysis of location plots from 1963 to 1977 
(Appendix A) shows this same general lack of data for these two regions. Statistical 
relevance must be considered when trying to calculate spatial scales with such 
sparsely located observations. 

Temporally uneven distributions are just as great. Over all, the data base 
ranges from a low of one profile for 1944 to a high of 1790 profiles in 1968. The 
number of yearly temperature profiles in the 15 year analysis period (1929 to 1991) 
is ranges from 592 to 925. When further subdivided into seasons, this greatly 
reduces the pool of temperature profiles for which statistical analysis may be 
conducted. Figure 4.3 and Appendix B provide examples of the temporal 
irregularities and sparseness of the observations. There are frequent 15 to 20 day 
gaps of ten observations or less per day. For an area as large as the Yellow Sea this 
does not provide an adequate data base to establish the inherent variability on a 


yearly basis. 
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Figure 4.2 (a) Winter Temperature Profile Locations from MOODS Data Base From 
1963 to 1977, (b) Locations of Temperature Profiles During Spring 1969, a Period 
Illustrating Maximum Data Density in a Given Year. 


Number of profiles 
N 
163 | 


? \ i | | Wt, | 
‘i Ae , A 
20 40 60 80 100 120 140 160 180 
Julian day 


Figure 4.3 Number of Temperature Observations During Winter and Spring 1977, 
Ilustrative of Typical Data Density in Any Given Half Year. 


18 


Spatial and temporal irregularities along with the shear lack of data must be 
carefully weighed in order to avoid statistically-induced variability. Originally, it 
was hoped that there would be sufficient data to compute horizontal autocorrelation 
scales on a yearly basis and determine a trend. After analyzing the location and 
histogram plots (Appendix A and B), it is obvious that statistically significant results 
can not be obtained on a yearly basis. Analysis was done on a mean seasonal basis 


using the data for all years. 
B. CLIMATOLOGY OF THE YELLOW SEA TEMPERATURE FIELD 


1. General 

Building a horizontal temperature field climatology from the MOODS data 
base is not a straightforward task and in shallow water the problems are magnified. 
Most of the problems relate to lack of data in contrast to temporal and spatial 
variability. Without enough data single observations force artificial gradients to 
appear. Seasonal time scales are thought to be more important than spatial scales but 
shallow water mitigates this by introducing much smaller temporal scales. In 
addition, the sloping bottom of the Yellow Sea poses difficulties with interpolating 
techniques. 


2. Method for Establishing Temperature Climatology 

Four main steps were used in establishing the temperature climatology: (1) 
seasonally binning the data, (2) interpolating the temperature profiles for the desired 
depth, (3) gridding the temperature with the Minimum Curvature Spline Interpolator 
(MISP) algorithm and (4) determining the desired grid spacing. The entire MOODS 
data base was binned by three month seasons as described earlier and then the 
temperature profiles were interpolated for the specified depths. Climatology 
building in deep water is distinctly different to that of shallow water. To compute 


deep water climatologies several submodels are constructed in the vertical (Figure 
4.4a). These submodels extend temperature profiles that do not extend to the sea 
floor (or to a given depth level) by interpolating values from profiles that do extend 
each submodel level. 

In shallow waters having sloping bottoms, if vertical submodels are used 
(Figure 4.4b), some areas (hashed area) are not covered by any submodels and 
therefore no climatology can be constructed for the area (Haeger, 1991). In order to 
compensate for this problem a vertical interpolating scheme that produces 
horizontal sloping surfaces was used (Figure 4.5). Horizontal temperature 
climatologies were constructed for the surface, 50% and 80% depth of the water 
column. Mid-water properties were expected to be modeled by the 50% depth 
whereas, the 80% depth was chosen since it was near the bottom but provided 
enough height for small scale depth variations. The surface temperatures in this case 
were actually any observations from the surface to 2 m below. 

After the temperature profiles were interpolated for the specified depths the 
temperature values were then horizontally gridded using the Minimum Curvature 
Spline Interpolator (MISP) algorithm. The MISP can be broken into three stages: 
partial grid computation, regional grid computation and final computation. The 
partial grid computation consists of binning the input data into bins equal to the 
specified output grid size. All data falling within a distance of half a grid cell of a 


grid point (i, 7) were averaged to form an averaged Xp y (location) and Tr, 
(temperature) set associated with that grid cell. At this stage T, is not situated at the 
gridpoint (i, j) but is located at the position (Xj, yi) . After all data have been read, 
a ‘partial grid’ has been formed and consists of grid points with and without 


associated (T..x.. y..) values. 
lady ay 
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The second stage of the gridding process consists of computing a regional 


grid. The (T,, x,» y,;) sets obtained for a grid cell in the first stage are used as input. 


(a) (b) 
top model 





top model 


middle model 


Figure 4.4 The Geometries of Vertical Submodel in the Construction of Climatology 
Temperature Profiles in the (a) Deep Water and (b) Continental Shelf. 


surface 
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Figure 4.5 Depiction of the Depth Levels that Horizontal Temperature 
Climatologies were Constructed for in the Yellow Sea. 


These values are first smoothed by rebinning to a grid three times the size of the 


original grid. A grid cell may or may not have a value at this stage. The minimum 
curvature spline interpolation routine was used to interpolate temperature (7;,) 


values to the gridpoints (i, j) centered at the grid cell resulting in the completion of 
the regional grid. 
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The third stage is to adjust the mean temperature (T;,) values of the partial 
grid to their respective gridpoints (i, j). First, the regional grid is interpolated back 
to the original grid size with a bi-cubic spline to match the partial grid. Then the 
gradient computed for each regional grid cell to the associated partial grid cell. The 
new partial grid is then merged with the regional grid. The merging consists of using 
the partial adjusted grid values where they exist at a gridpoint and the regional grid 
values where the partial grid cells are empty. This allows the input data to dominate 
in data rich areas and the regional data to dominate in data sparse areas (Haeger, 
personal communication). 

Sometimes a climatology field is sensitive to the selection of the grid size. 
Too small a grid spacing produces false gradients appearing as bull’s-eyes. Yet, if 
the spacing is too large, features that do exist may not appear. Grid spacings of 20 
minutes and 30 minutes were analyzed in this study. As expected, the contours 
associated with a 30 minute grid are smoother (Figure 4.6a) than the 20 minute grid 
(Figure 4.6b) causing mesoscale features to be lost. In general, the 20 minute grid 
better depicts the findings of Elms (1989) for the Yellow Sea and hence this grid will 
be used for the rest of the study. 


3. Yellow Sea Mean Temperature Fields 

The surface temperatures for all four seasons (Appendix C) show that 
temperatures increase toward the south. During winter, freezing temperatures occur 
in the Bohai Sea region which result in surface ice. In the Yellow Sea winter 
temperatures are significantly higher, with readings in the high teens (14-18 °C) near 
the shelf break. As the year progresses towards summer, the horizontal temperature 
gradient decreases. Spring temperatures range from 6 °C in the Bohai Sea, 14 °C in 
the central Yellow Sea and 24 °C near the Kuroshio Current. By the ttme summer 
arrives the SST of the entire Yellow Sea is above 20 °C and is almost uniformly the 


Jags 


same. With the onset of cooler fall conditions the north-south gradient reappears. 
The Bohai Sea shows the presence of a pool of warm (18 °C) water while the Yellow 
Sea begins to establish its winter gradient with lower temperatures (14 °C) to the 


north and a steady increase to the shelf break (24 °C). 
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~ 4.6 Climatological SST Field on (a) a 30 Minute Grid and (b) a 20 Minute 
rid. 

The horizontal temperature contours at 50% and 80% water column depths 
tend to mimic each other from season to season. The winter subsurface 
climatologies closely mirror that of the surface. The largest difference is in the shelf 
break area where surface temperatures are slightly higher. On the other extreme, 
summer shows the most contrast between surface and subsurface temperatures. The 
surface temperature field shows almost uniform values throughout whereas the 50% 
and 80% depth contours show a cold core in the middle of the Yellow Sea. Both 50% 
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and 80% depth levels show this cold core beginning to form in spring and 
dissipating in fall. 
C. SPACE / TIME SORTING 

At present, the spatial and temporal scales for the OTIS model are artificially 
prescribed as 100-250 km (Phoebus, 1988) and 20 days (Clancy et al., 1989) for the 
region of Japan and the South China Sea. Since synoptic temporal and horizontal 
scales are usually smaller than 30 days and 300 km (Mooers and Robinson, 1981), 
it was decided that a matrix with a maximum temporal lag of 30 days with a one day 
increment and the maximum spatial lag of 300 km with a 10 km increment would 
be used. Each matrix bin is associated with a specific temporal and spatial lag. Other 
matrix sizes were tested. It was found that smaller temporal and spatial lags 
excluded too much data while larger lags did not necessarily improve the results. 

As mentioned earlier, an important feature of the Yellow Sea thermal field is 
the existence of the CWYS (Li and Yuan, 1992). This water mass is located in the 
central bottom of the Yellow Sea. In the winter it merges with the surface water but 
from April to November it is a distinct water mass with temperatures of 4-5°C. 
During the summer the overlying waters are 10-15°C warmer. The isotherms are 
closed curves (Figure 4.7), suggesting the basin may be treated as an isolated system 
(Guan, 1994). Therefore, it is reasonable to assume that the spatial autocorrelations 
of temperatures inside the Yellow Sea basin are isotropic, 1.e.,the autocorrelations 
are independent of direction (N-S or E-W) and depend only on the distance between 
two temperature locations. This assumption of horizontal isotropicity greatly 
simplifies the computation. 

After the matrix size was defined and the isotropic assumption made, 


anomaly pairs were placed within the temporal and spatial bins. To accomplish this, 


it is necessary to take each observation,7, find the closest gridpoint and compute 
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the anomaly, 7,’, by subtracting the climatological value at this gridpoint ie from 
the observation, i.e., 7° = T,- The For every individual temperature anomaly, T,’, 


it is paired with the other data points, 7’, within the data base. The temporal and 


38° 


36° 





Figure 4.7 Temperature Distribution in Summer (April-November) of the ‘Cold 
Water Mass of the Yellow Sea’ (CWYS) Near the Sea Floor (from Su and Weng, 
1994). Isotherms in Degrees Celsius. 


spatial difference or lag is then calculated between the two anomalies. The anomaly 


pairs (T\’,7,’) are then placed in the corresponding temporal and spatial lag bin. 


0’ ° Oo 
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Three separate data runs were conducted. The first data run was to use 
separate seasons and individual years from 1963 to 1977. The second run was to use 
all the MOODS data binned by season. There was no delineation between years; 
therefore an observation from 23 January 1963 could be paired with an observation 
from 24 January 1969 and only have a lag time of one day. The third run also 
included the entire data set binned seasonally but each observation retained its 
year.Hence, it was paired with observations of the same year, then sorted into the 
matrix. For simplicity the different input types will be referred to as Type I, Type IJ, 
and Type LI. 

The number of observations 1n the data set (V) will determine the number of 
pairs that can be binned. The number of pairs (S) is computed by 
_ N(N-1) 

Z 
In Type I and Type II data the computation S is straight forward but for Type II the 


S (4.1) 


value of S is the summation of S for each individual year. As to be expected, Type I 
has the fewest pairs in each bin. It ranges from one pair to a high of 2000 pairs. In 
contrast Type I ranges from approximately 750 to over 24000 pairs and Type Hil 
from approximately 40 to the 4000. Figure 4.8 shows the pair distribution of Type I 
and Type I. It can be easily seen that the Type I scheme (Figure 4.8a) produced 
three main (green) concentrations of profile pairs, the other areas lack substantial 
data. Type II (Figure 4.8b) has a much higher and more even distribution. From this 
lack of data in Type I it was determined that statistically relevant information could 
not be obtained. This study will analyze the advantages and disadvantages of Type 
II and Type IL data. 
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Figure 4.8 Contour Plots of Time-Distance Lag Pair Distribution (a) TYPE I for 


Spring of 1969 and (b) TYPE II for Spring. 
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D. AUTOCORRELATION FUNCTION 

After the anomaly pairs have been spatially and temporal sorted, an 
autocorrelation function (nj (m,n)) for each bin with temporal lag nAr and 
horizontal lag mAy (At = 1 day, Ay = 10 km) identified as the bin (m, n) is computed 
by 


oe, ts 


bin (m, n) (4.2) 


Nn (m,n) = ; ; 
| y To | d to 
bin (m, n) bin (m, n) 


Figure 4.9 shows the temporal variation of autocorrelation ( (m,n) 7 at the 
surface for the winter season of Type II and Type III data. The three curves on each 
graph show the dependency of n(m,n) on temporal lag nAt for three different 
spatial lags: m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and m = 15 (150 km 
lag, dots). Such a decrease of n (m,n) with nAt reveals approximately a Gaussian 
type relation. On both graphs there is no unity or a correlation of one at zero lag. This 
occurs since the covariance estimates at zero temporal lag did not include products 
of each data point with itself and also because substantial subgrid noise was present 
(White et al., 1982). 1 (0,0) is greater for Type II than the Type II computations 
which is expected since temperatures of the same year are being compared with each 
other. Type I plots (Figure 4.9a) are much smoother than the Type IH plots (Figure 
4.9b). This high frequency temporal variability may be induced by a smaller data 
base and by forcing mechanisms like meteorological, tidal and freshwater runoff 
that were averaged out in the Type II computations but are clearly recorded in the 
Type If calculations. It is unclear how many or what type of natural forcing 


mechanisms are creating this variability or if it is an artifact of having too little data 
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that artificially induces variability. In each case the e-folding time for zero lag 
distance is approximately 13 days. 

The spatial autocorrelations showed similar characteristics (Figure 4.10). 
Type III had higher autocorrelations at zero lag. Variability is large in the Type Il 
scheme and the reason for this is still unclear. The e-folding distance for Type II data 
was roughly estimated from Figure 4.10 to be approximately 125 km while for the 
Type III it is 150 km. Since the e-folding scales are similar for both computations 
(Type II and Type III) and Type II exhibits statistical stability, the Type Il 
computations were considered the more appropriate choice in calculating the 


decorrelation scales. 


E. E-FOLDING DECORRELATION SCALES OF THE YELLOW SEA 
THERMAL FIELD 


With regard to the surface and subsurface thermal indicator parameters, the 
e-folding temporal and spatial decorrelation scales can be obtained from the 
autocorrelation coefficient function, nN (m,n) . The variations of n (m,n) verses the 


horizontal spatial and temporal lags (mAy and nAt) reveals a Gaussian-type 


dependency (e.g., Figure 4.9). The autocorrelation matrices n () (m,n) are fitted by 


the Gaussian function 


—A*(mAY)?— Cy (nAt)’ 
n (m,n) = oe 0 ¢ (4.3) 


where k denotes the different levels (k =1 for the surface, k =2 for the 50% water 
column depth, and k =3 for the 80% water column depth). The regression 


coefficients, A;' and C;', are the e-folding horizontal spatial and temporal 


decorrelation scales. For a given time lag nAt, the coefficient Am is interpreted as 


an e-folding distance, m,Ay, for the level k 
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Figure 4.9 Temporal Variations of Autocorrelation Function 7 (m,n) for (a) Type 
If and (b) Type If Data at the Surface.The Three Curves on Each Graph are the 
Dependency of n (m,n) on Time Lag nAt for m = 0 (no lag, dash-dot), m= 1(0km 
Jag, solid), and m = 15 (150 km lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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Figure 4.10 Spatial Variations of Autocorrelation Function n (m,n) for (a) Type Il 
and (b) Type III at the Surface. The Three Curves on Each Graph are the 
Dependency of n (m, n) on Spatial Lag (mAy) for n =0 (no temporal lag, dash-dot), 
n= 1 (1 day lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the 
MOODS Temperature Data Base for the Years 1929 to 1991. 
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n i (m,n) 


-A,”(mAy) = In = -] (4.4) 


k 
n (0,n) 
where the autocorrelation function reduces from zero spatial lag, 7” (0,n), to 
m.Ay spatial lag, 7” (m,n) = ely (O, 7). 
Similarly, for a given spatial lag (mAv¥) , the coefficient C, ' is interpreted as 
an e-folding temporal separation (n,At) for the level k 
k 
—C, (n.At)~ = In————— = - 1 (4.5) 
: (k) 
n> ~ (m, 0) 
where the autocorrelation function reduces from M7 (4) (m,0) to 


n° (m,n) = en (m,0). 


F. SIGNAL-TO-NOISE RATIO 


In applying the optimum interpolation for the OTIS model or an optimal 
network design, certain statistics of the field of interest such as the variance of both 
the signal and the noise must be known. The signal variance for the bin (m, n) 1s 


estimated by 


2 _ Ad 
Ea ee ie, 
bin (m, n) 


(4.4) 


The noise variance is derived from two sources (Sprintall and Meyers, 1991), 
geophysical error and instrument error. Geophysical error is unresolved thermal 
variability with scales smaller than the typical time and distance scales between two 
temperature profiles. Instrumental error can be shown to be negligible compared to 
geophysical noise (White and Bernstein, 1979). These noise sources inhibit the 


o2 





detection of large-scale variability and can be an important factor in establishing the 
sampling density necessary to detect the signal (Sprintall and Meyers, 1991). 

In calculating the autocorrelation value of the first bin (zero to 1 day and zero 
to 10 km lag) the temperature anomalies were never paired with themselves. The 


values of the autocorrelation coefficient function at zero temporal and spatial lag 


[n (0, 0) ] were extrapolated from the first bin into the origin. Therefore, n (0, 0) is 
not equal to 1.0. The relationship between the temperature signal and the noise 
signal can be described by the signal-to-noise parameter (Sprintall and 


Meyers, 1991) 
0) 
es A ae ee) (4.5) 
1-7n(0,0) oOo 


n 
where o is the noise variance and Go is the signal variance. A represents the signal- 
to-noise ratio or level. The larger the A, the less noise that exists in the bin. 
Theoretically, if there is no noise, n (0,0) = 1 and o = Q. Therefore 4 will be ~. 


If > 2, the ratio of the signal variance, O°, to the noise variance, a is greater than 


4 which is considered good by White et al. (1982) and Sprintall and Meyers (1991). 
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V. DISCUSSION 


A. TEMPORAL AND SPATIAL DECORRELATION SCALES 

All temperature profiles from 1929 to 1991 were used in the calculation of 
temporal and spatial decorrelation scales. Delineation of separate years within the 
autocorrelation computations was not permitted (Type II averaging) which averaged 
the various forcing mechanisms prevalent in this shallow water basin. The computed 
autocorrelations for each season were then fitted to a Gaussian function of the form 


of (4.3). The regression method was used to calculate the spatial and temporal 


; : =I = ; : 
decorrelation scales at three different levels, 4, and C,’, respectively. Since an 
isotropic assumption concerning the spatial temperature distribution was made, 
there was no directional dependency. This means that decorrelation scales were 


assumed to be the same throughout the Yellow Sea basin. Table 5.1 presents these 


decorrelation scales and signal-to-noise ratios for winter and summer. 


| Season Depth Temporal Spatial Signal-to 
days km noise ratio(A 


Summer 





Table 5.1: Seasonal Temporal and Spatial Scales 


The signal-to-noise ratio for both winter and summer are higher than the 
signal-to-noise ratio for the deep waters of the Eastern Equatorial Pacific (~ 1.0) as 
reported by Sprintall and Meyers (1991). This infers that the temperature signal in 
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this study is stronger (more stable) than the background noise. In both seasons the 
smallest signal occurs at the 50% water column depth. This is expected since this is 
the transitional (thermocline) depth between the surface layer and the deep layer, 
where meteorological and topographic effects occur. The surface signal-to-noise 
ratio was smaller in winter than in summer. 

Winter temporal and horizontal decorrelation scales are the smaller than 
Summer and exhibit nearly uniform values at all three depths. The temporal scales 
vary only by a half a day and the horizontal scales by only 14 km. This uniformity 
is expected considering the oceanographic conditions of the region. Sea surface 
icing occurs to the north in the Bohai Sea and the northern Yellow Sea and along the 
northern coastal regions of Korea and China. Temperatures range from freezing in 
the north to 20°C in the south along the shelf break. This large temperature range 
will inherently make horizontal decorrelation scales small. Vertically, winter 
experiences the greatest mixing (Nester,1977) with monsoonal winds at their 
strongest throughout the entire season. With a large part of the Yellow Sea having 
depths less than 50 m, mixing reaches to the bottom creating isothermal profiles. 
Thus, decorrelation scales will be similar from the surface to the 80% water column 
depth. 

Summer surface spatial decorrelation scales are 100 km longer than winter. 
This is reasonable since this region experiences large solar heating in the summer. 
The summer climatology (Appendix C, Figure C.1) shows the SST to be nearly 
constant throughout the region. The surface temporal scale is less by five days than 
at deeper depths. This may be due to a combination of factors. Besides the large 
solar heating that occurs in the summer, the summer monsoonal winds are much 
weaker than the winter winds. Therefore, vertical mixing is not as great as in winter. 
A strong thermocline is present in summer. Thus, the upper layer water is affected 


by the meteorological synoptic events rather than the entire water column as is in 
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winter. The temporal scale at the surface should be shorter in summer than in winter 
(Table 5.1). In the deeper layers (50% and 80% depths) temperatures are not directly 
affected by the meteorological events of the surface and experience longer temporal 
scales. The 80% water column depth has the longest temporal scale which could be 
directly related to the ‘Cold Water Mass of the Yellow Sea’ (CWYS) (Li and Yuan, 
1992) that exists throughout the summer in the middle of the Yellow Sea. Spatial 
scales are close to 100 km smaller than that at the surface. Without the extreme 
vertically mixing, the CWYS retains some its winter water mass characteristics from 
which it is derived (Liu et al., 1992) and remains a distinct cold water mass that 


resists mixing. This will cause smaller decorrelation scales. 


B. OBSERVATIONAL NETWORK DESIGN 

The temporal and horizontal decorrelation scales provide useful information 
for an optimal design of an observational network. The principle is to minimize both 
instrumental and geophysical errors. Since noise due to instrumental error in XBT 
measurements is usually only about 0.1°C (Barnett and Bernstein, 1980), and even 
smaller in CTD measurements, instrumental error can be neglected. This implies 
that the reduction of noise must be accomplished by the reduction of geophysical 
error. This is usually done by increasing the sample density. Sprintall and Meyers 
(1991) suggested that an optimal sampling density of two or three samples per 
decorrelation scale would be adequate to resolve thermal variability. 

This would mean that spatially, any future temperature measurement in the 
winter may be conducted at 50-80 km intervals with the knowledge that the 
subsurface features will also be adequately sampled. The summer season will 
require more consideration. If only surface features are considered, a larger 


measurement interval of 80-125 km may be used. In order to provide full coverage 
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of subsurface features a 50-80 km interval is recommended.Temporally, a 4-6 day 


sampling interval will provide the minimal coverage for both seasons. 


C. UNCERTAINTIES IN THE TRANSITIONAL SEASONS 

The transitional seasons, spring and fall, showed unusual temporal and 
spatial decorrelation scales (Table 5.2). The scales were significantly larger than in 
winter or summer. The spatial surface scales computed are 521 km in the spring and 
757 km in the fall. It might be expected that the transitional surface scales be longer 
than winter and summer due to the lack of meteorological synoptic forcing or the 
tendency for the temperature contours (Appendix C), especially in spring and to a 
lesser degree in fall, to rap around or follow the coastline around the Yellow Sea 


basin. This could mean that an eastern coastal temperature profile would be more 


correlated to a western coastal profile 300 km away than it is to the profiles in the 


Depth Temporal Spatial Signal-to 
days cr noise ratio(A 


Spring 





Table 5.2: Seasonal Temporal and Spatial Scales 


middle of the basin that are closer. In either case, the extremely large scales are 
unrealistic if not for the only reason that the enclosed Yellow Sea basin itself is 
approximately 450 km in the east-west direction and over 700 km in the north-south 


direction. 
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The temporal scales are approximately six days longer than in winter and 
summer. There was no evidence to support these larger values. Since the onset or 
cessation of seasonal changes vary from year to year and this study specifically 
defined April, May and June as spring and October, November, and December as 
fall, statistical errors are expected. It is highly recommended that future studies 
examine different month groupings for these transition seasons. 

Vertically both the temporal and spatial scales decreased with depth. The 
larger change is in the spatial scales. In spring there was a 300 km and in fall a 500 
km difference from the 50% and 80% depths to that at the surface. As stated earlier, 
these findings seem highly improbably and the irregularities of the MOODS data 


and the determination of season must be examined further. 


D. LIMITATIONS 

The purpose of this study was to calculate the seasonal temporal and spatial 
decorrelation scales to be used by OTIS. Scales were obtained but the validity of 
these findings must be examined and a potential user discouraged from blindly 
implementing these scales into various models. 

The isotropic assumption of uniform conditions throughout the basin was 
based on the existence of closed isotherms for the majority of the time in the Yellow 
Sea basin. This made computations feasible and provided a general reference for 
temporal and horizontal scales. In reality, there probably exists a similar latitudinal 
and longitudinal dependence as in deep water (White et al., 1982; Sprintall and 
Meyers, 1991). This is due to the fact that isotherms are not closed year round, the 
basin 1s elongated in the north-south direction and coastal regions typically have 
different cross-shore and along-shore scales. Thus, future studies should examine 


the possibilities of calculating two directional scales. 
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The convention used in defining seasons can have a large statistical impact. 
The convention used in this study was recommended by the Naval Oceanography 
Office. It was found that spring and fall had unrealistic surface scales which may be 
a direct result of the seasonal binning. A more realistic approach may be to define 
spring and fall as two months or less in duration as opposed to three used in this 
Study. 

The MOODS data base is a limiting factor of itself. In this study, 
decorrelation scales that pertain to the entire Yellow Sea basin were desired. Yet, the 
data were biased in location (Appendix B). Most of the profiles were located in the 
eastern region of the sea, with significantly fewer in the Bohai Sea and along the 
coastline of China. Temporally the profiles were irregular (Appendix A) from 
month to month and year to year. This irregularity led to the analysis of three 
different computational schemes. 

The first scheme, Type I, calculated autocorrelations on a yearly basis. 
Values were erroneous since not enough data existed to produce statistically 
relevant results. Type II calculations used the entire MOODS data base but no 
delineation was permitted between years. This scheme averaged the high frequency 
variability such as meteorological events, tidal forcing and river runoff. This scheme 
showed a definite seasonality (Appendix D) both temporally and spatially. Type Il 
calculations also used the entire MOODS data but this time temperatures were 
paired with only temperatures of the same year. Although the autocorrelations 
computed from Type II and Type III schemes are similar, the Type Il scheme 
produced large variability both temporally and spatially (Appendix D). It is not 
known if the variability is statistically produced by a smaller data base or are the 
result of the many shallow water forcing mechanisms. Therefore, more effort is 


needed to determine the best use for both schemes. 
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The last limitation to be discussed is the fitting function. Even though it is 
customary to fit the calculated autocorrelations to a Gaussian curve, as in the OTIS 
model, it may not be the most accurate procedure for shallow water. The winter and 
Summer autocorrelations rapidly decreased in time and space. The spring and 
summer curves did not. In future studies it may be better to use a combination of a 


linear function and a Gaussian curve. 
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VI. CONCLUSION 


The goal of this study was to compute seasonal temporal and spatial 
temperature autocorrelation scales for the Yellow Sea. These scales are an integral 
part of the accurate running of the OTIS model which at this time is designed for 
deep water thermal analysis. Decorrelation scales were calculated for three depths, 
i.e., the surface and at 50% and 80% of the water column depth, and for four seasons. 
Variations were found to be both seasonally and depth dependent. 

The temperature data were comprised of approximately 36,000 observed 
profiles from the MOODS worldwide data base. All profiles were located on the 
continental shelf and spanned the time period of 1929 to 1991. The widely irregular 
distribution of temperature profiles in time and space made yearly analysis 
Statistically irrelevant. Two time-organized schemes were analyzed: (1) 
temperatures were paired with no delineation between years, referred to as Type I 
data and, (2) temperatures paired with temperatures of the same year, referred to as 
Type III data. The Type II scheme was used in the calculation of the temporal and 
spatial scales since it produced a larger data set and averaged the high frequency 
variability due to individual forcing mechanisms. 

The winter season had the smallest temporal and spatial scales of the four 
seasons. This is attributed to the strong cold monsoonal winds which produced the 
strongest vertical mixing of all the seasons and the extreme temperature range 
extending from freezing in the north to 20°C in the south along the continental shelf. 
The temporal decorrelation scales approximated to 15 days and spatial decorrelation 
scales to 165 km at all depths. 

The summer season had the next smallest scales, but unlike winter it showed 


some vertical variability. The scales range from 12.3 days and 251 km at the surface 
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to 17.2 days and 157 km at 80% depth. The approximately 100 km difference 
between the levels is due to the lack of vertical mixing and the presence of CWYS. 

An optimal network design was suggested. It was found that a 50-80 km 
interval will provide adequate sample density to discern winter features. The 
Summer season exhibits a surface and subsurface minimum density interval. To 
discern only surface features a 80-125 km interval is adequate but for complete 
depth coverage a 50-80 km interval should be used. Temporally, a 4-6 day sampling 
interval will provide coverage for both seasons. 

The use of these scales in the OTIS model must be cautiously considered. The 
inherent problems in the MOODS data base, irregularity in space and time, the 
organization of the data base into seasons and the neglect of different forcing 
mechanisms all have an effect on the calculation of the temporal and spatial scales. 
More temperature profiles and analysis is required in order to better refine these 


scales and improve the OTIS model. 
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APPENDIX A SEASONAL LOCATION PLOTS OF MOODS 
TEMPERATURE PROFILES FROM 1963 TO 1977 
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Figure A.2 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 


for 1964. 
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Figure A.4 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 


for 1966. 
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Figure A.5 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1967. 
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Figure A.6 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1968. 
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Figure A.7 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 


for 1969. 


Winter Spring 






: ecocccascbeos 


CF © oe & AAPoR oot qr 
i 
i 


Latitude (N) 


Coes ceeecoqas 


e 
ep, *? 
o°8 
* 
e 

oy e@ 
@ 





118 120 122 124 126 128 118 120 122 124 126 128 


Summer Fall 


eeqececotoe 
5 a he 33° 


Latitude (N) 





118 120 122 124 126 128 118 120 122 124 126 128 
Longitude (E) Longitude (E) 


Figure A.8 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1970. 
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Figure A.9 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 


for 1971. 
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Figure A.10 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 


for 1972. 
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ligure A.11 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1973. 
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Figure A.12 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1974. 


60 


Winter 


Ww 
oN 


eeeeeseceange 


eo ¢@ 
eeeee Pee 


Latitude (N) 
Nh Nh w 7) 
o>) @ © Nh 


Nh 
BS 


118 120 122 124 126 128 


Summer 


TMs 


eccccceevee 


#0) 
BAS 


o 
@eeseeeeecoeee 
e 


Latitude (N) 
w 
NO 


Ww 
© 


© g ere cgte 
ee *e 


NO 
co 


nN 
oO 


Nh 
oS 


118 120 122 124 126 128 
Longitude (E) 


Figure A.13 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 


for 1975. 
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Figure A.14 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 


for 1976. 


62 


Winter Spring 


34 eeee ccc, = 4 
meee, $3 


ies ial 


e¢ 
®e% 


@eee0et? #0 
Vd 


Latitude (N) 
Gd 
NO 





eheecccac ee of: Co 
30 Steger cc's 5 
28 H 
: 
26 ° 
24 
118 120 122 124 126 128 118 120 122 124 126 128 
Summer Fall 


Bp Povecde 


Latitude (N) 
ie) 
i) 





30 or fcce eee se 
are 
28 
26 
24 
118 120 122 124 126 128 118 120 122 124 126 128 
Longitude (E) Longitude (E) 


pene A.15 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
or 1977. 
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APPENDIX B TEMPORAL DISTRIBUTION OF MOODS 
TEMPERATURE PROFILES FROM 1963 TO 1977 
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Figure B.1 The Number of MOODS Temperature Profiles During 1963 in the 
Yellow Sea. 
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Figure B.2 The Number of MOODS Temperature Profiles During 1964 in the 
Yellow Sea. 
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Figure B.3 The Number of MOODS Temperature Profiles During 1965 in the 
Yellow Sea. 
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Figure B.4 The Number of MOODS Temperature Profiles During 1966 in the 
Yellow Sea. 
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Figure B.5 The Number of MOODS Temperature Profiles During 1967 in the 
Yellow Sea. 
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Figure B.6 The Number of MOODS Temperature Profiles During 1968 in the 
Yellow Sea. 
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Figure B.7 The Number of MOODS Temperature Profiles During 1969 in the 
Yellow Sea. 
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Figure B.8 The Number of MOODS Temperature Profiles During 1970 in the 


Yellow Sea. 
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Figure B.9 The Number of MOODS Temperature Profiles During 1971 in the 
Yellow Sea. 
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Figure B.10 The Number of MOODS Temperature Profiles During 1972 in the 
Yellow Sea. 
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Figure B.11 The Number of MOODS Temperature Profiles During 1973 in the 
Yellow Sea. 
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Figure B.12 The Number of MOODS Temperature Profiles During 1974 in the 
Yellow Sea. 
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Figure B.13 The Number of MOODS Temperature Profiles During 1975 in the 
Yellow Sea. 
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Figure B.14 The Number of MOODS Temperature Profiles During 1976 in the 
Yellow Sea. 
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Figure B.15 The Number of MOODS Temperature Profiles During 1977 in the 
Yellow Sea. 
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APPENDIX C YELLOW SEA SEASONAL MEAN TEMPERATURE 
FIELDS 
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Figure C.1 Yellow Sea Seasonal Mean Temperature Field for the Surface on a 20 


Minute Grid. Data were Obtained from the MOODS Temperature Data Base for the 
Years 1929 to 1991. 
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APPENDIX D AUTOCORRELATION FUNCTIONS 
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Figure D.1 Temporal Variations of Autocorrelation Function n (m,n) for Type Ul 
Data at the Surface.The Three Curves on Each Graph are the Dependency of 

1 (m, n) on Time Lag nAt for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m= 15 (150 kmlag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.2 Temporal Variations of Autocorrelation Function n (m,n) for Type I 
Data at 50% Water Depth. The Three Curves on Each Graph are the Dependency of 
1) (m,n) on Time Lag nAt for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.3 Temporal Variations of Autocorrelation Function n (m,n) for Type I 
Data at 80% Water Depth. The Three Curves on Each Graph are the Dependency of 
N (m, n) on Time Lag nAt for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m= 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.4 Temporal Variations of Autocorrelation Function n (m,n) for Type Il 
Data at the Surface. The Three Curves on Each Graph are the Dependency of 

N (m, n) on Time Lag nAt for m =0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m= 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.5 Temporal Variations of Autocorrelation Function n (m,n) for Type IL 
Data at 50% Water Depth. The Three Curves on Each Graph are the Dependency of 
1 (m,n) on Time Lag nAt for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m= 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.6 Temporal Variations of Autocorrelation Function n (m,n) for Type Ii 
Data at 80% Water Depth. The Three Curves on Each Graph are the Dependency of 
1] (m,n) on Time Lag nAt for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.7 Spatial Variations of Autocorrelation Function n (m,n) for Type II at 
the Surface. The Three Curves on Each Graph are the Dependency of 7 (m,n) on 
Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day lag, solid), 
and n = 15 (15 day lag, dots). Data were Obtained from the MOODS Temperature 
Data Base for the Years 1929 to 1991. 
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Figure D.8 Spatial Variations of Autocorrelation Function n (m,n) for Type I at 
50% Water Depth. The Three Curves on Each Graph are the Dependency of 

1 (m,n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day 
lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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Figure D.9 Spatial Variations of Autocorrelation Function n (m,n) for Type II at 
80% Water Depth. The Three Curves on Each Graph are the Dependency of 

1 (m,n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day 
lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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Figure D.10 Spatial Variations of Autocorrelation Function n (m,n) for Type II at 
the Surface. The Three Curves on Each Graph are the Dependency of n (m,n) on 
Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day lag, solid), 
and n = 15 (15 day lag, dots). Data were Obtained from the MOODS Temperature 
Data Base for the Years 1929 to 1991. 


94 


Winter Spring 


Autocorelation coefficient 





0 100 200 300 


Summer Fall 


Autocorrelation coefficient 





0 100 200 300 0 100 200 300 
Distance lag (km) Distance lag (km) 


Figure D.11 Spatial Variations of Autocorrelation Function n (m,n) for Type II at 
50% Water Depth. The Three Curves on Each Graph are the Dependency of 

1 (m,n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day 
lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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Figure D.12 Spatial Variations of Autocorrelation Function n (m,n) for Type Il at 
80% Water Depth. The Three Curves on Each Graph are the Dependency of 

1 (m,n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day 
lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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